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METHODS AND SYSTEMS FOR 
TIMING-DRIVEN CIRCUIT SYNTHESIS 



5 Background of the Invention 

Field of the Invention 

Aspects of the present invention generally relate to computer aided engineering 
of logic circuits. More particularly, embodiments of the present invention relate to 
timing optimization of logic circuits. 

10 

Description of the Related Technology 

In general, logic optimization is classified into two categories, two-level logic 
optimization and multi-level logic optimization. 

Two-level optimization deals with the optimization of combinational logic 
15 circuits, modeled by two-level "sum of products" expression forms, or equivalently by 

tabular forms such as implicant tables. Two-level logic optimization has a direct impact 
on programmable logic arrays (PLAs) and macro-cell based programmable logic 
devices (CPLDs). 

Combinational logic circuits are very often implemented as multi-level networks 
20 of logic gates. The fine granularity of multi-level networks provides several degrees of 
freedom in logic design that may be exploited in optimizing area and delay as well as in 
satisfying specific constraints, such as different timing requirements on different 
input/output paths. Thus, multi-level networks are very often preferred to two-level 
logic implementations such as PLAs. The unfortunate drawback of the flexibility in 
25 implementing combinational functions as multi-level networks is the difficulty of 

modeling and optimizing the networks themselves. The need of practical synthesis and 
optimization algorithms for multi-level circuits has made this topic of high importance 
in VLSI CAD. 

Multi-level logic optimization is fi*equently partitioned into two steps. In the 
30 first step, a logic network is optimized while neglecting the implementation constraints 
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on the logic gates and assuming rough models for their area and performance. This 
procedure is usually referred to as technology independent logic optimization. In the 
second step, one takes into consideration the constraints on the available gates (e.g., K- 
LUTs in FPGAs) as well as the detailed area and delay models of these gates. This step 
5 is the so-called technology dependent logic optimization or technology mapping. The 

discussion hereinbelow addresses the technology independent logic optimization 
problem, and, in particular, the timing-driven logic resynthesis problem. 

Several common operations that are used during the area-oriented multi-level 
optimization are as follows: 
10 L Common sub-expression extraction 

By extracting common sub-expressions from a number of functions, the circuit 
area is reduced. However, the better the area saving, the more places the sub-expression 
fans out to, which could degrade the circuit performance. 

2. Resubstitution 

15 Resubstitution is similar to common sub-expression extraction and involves 

expressing a node in terms of another, if possible. 

3. Elimination 

Elimination involves removing, from the multi-level network, all occurrences of 
variables that represent the nodes which are eliminated. When all the internal nodes are 
20 eliminated, the operation is called collapsing. 

4. Decomposition 

The decomposition of an internal node function in a multi-level network 
replaces the node by two (or more) nodes that form a subnetwork equivalent to the 
original node. Decomposition is often performed on a node to split a complex function 
25 into two (or more) simpler functions. Small-sized expressions are more likely to be 

divisors of other expressions and may enhance the ability of the resubstitution algorithm 
to reduce the size of the network. 

5. Simplification using don't care conditions 

Simplification is used to find a compact representation for the Boolean function 
30 at every node. By removing the redundancies from a representation of a function, both 

the size and the depth can be reduced. In a multi-level network, the simplification at a 
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node needs to consider the structure of the logic around it. This gives rise to don't care 
conditions that can be exploited during node simplification. 

From the description of these operations, one can see the complex interaction 
between the circuit area and delay. In addition, the delay impact of a particular 
5 transformation applied on the same network often depends on the delay data (the arrival 

and required times). Since the delay data is imprecise at the technology independent 
stage, it is difficult to adapt the strategies used for area optimization to address the 
performance optimization issue. Because of this difficulty, many of the techniques 
developed to reduce the circuit delay use local transformations to make incremental 

10 changes to the logic. 

Timing optimization will now be discussed. One significant issue in 
restructuring a circuit is determining circuit regions that should be transformed. The 
most critical outputs and their transitive fanins are a natural choice. However, one 
problem with this approach is that after the most critical outputs have been optimized, 

15 outputs that were close to being critical before could become critical after optimization 

of the original critical paths. Moreover, optimizing only the most critical outputs by 
more than the needed arhount can also result in an unnecessary area penalty. Thus, 
some techniques optimize close-to-critical nodes along with the most critical nodes. 

Several conventional algorithms use an iterative refinement-based approach, 

20 where, in each iteration, a set of critical paths is identified and then the delays of a set of 
nodes are reduced so that the overall circuit performance is improved. These algorithms 
are differentiated in (i) how to determine in each iteration the set of nodes to apply the 
local transformation for delay reduction and (ii) the local transformation method itself 
Another conventional attempt at timing optimization takes a different approach 

25 based on clustering, partial collapsing and subsequent timing optimization. This 

approach is based on the premise that at a technology-independent level, in the absence 
of the target technology information and wiring delays, any delay model is inaccurate. 
Therefore, it assigns a zero delay to all the gates, thus treating all the input-to-output 
paths uniformly. However, whenever a signal crosses cluster boundaries, a delay of one 

30 unit is incurred. 
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Another existing approach first performs area optimization on a circuit to 
achieve to reduce the size of the circuit layout, and then incremental changes are made 
to the circuit to reduce its delay. This approach is particularly useful for layout-driven 
logic resynthesis, v^herein the timing correction step is performed incrementally to 
5 ensure the convergence of the iteration betv^een the layout design and the circuit 

resynthesis. 

A significant aim of the restructuring approaches discussed above is to generate 
a good multi-level structure of the circuit that will subsequently be mapped into a small 
delay implementation. These conventional approaches generally use simple, weak 
10 models to predict the circuit delay. As a result, the savings observed at the technology 

independent stage may not be evident after technology mapping of the optimized 
circuit. 

To alleviate this problem, researchers have extended the basic ideas of the 
technology independent optimizations to work on mapped circuits. Heuristics have 
15 been used to address the optimization of mapped circuits while taking into account the 

characteristics of the cell library. 

The Timing-Driven Logic Optimization section discussion below describes the 
performance optimization at the technology independent level and how this 
optimization impacts the subsequent technology mapping and physical design. 
W 20 With the rapid scaling of transistor feature sizes, integrated circuit performance 

U is increasingly determined by interconnects instead of devices. Interconnect delays are 

even more significant in PLD designs due to the extensive use of programmable 
switches. As a resuU, the delay between two logic blocks is highly dependent on their 
placement on the chip and the routing resources used to connect them. PLDs, such as 
25 those fi*om Altera, include several types of interconnects, including local, row and 

column interconnects. Local interconnects refer to the connections between logic 
elements (LEs) in the same logic array block (LAB). Row intercormects refer to the 
connections between LEs in the same row, but in different LABs. Column 
interconnects refer to the connections between LEs in different rows. The delay 
30 attributed to interconnects can be many times that of the logic element delay. Given 

such a high variation of different types of interconnect delays, it would be almost 
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impossible to perform accurate timing optimization during synthesis without proper 
consideration of the layout result. That is why layout-driven synthesis is considered to 
be an important problem area in high-performance PLD designs. 

The layout-driven synthesis problem has proved to be difficult to solve due to 
the mutual dependency nature of the logic synthesis and layout design. In general, there 
are two approaches to integrate logic and layout synthesis. One approach is to employ a 
highly iterative design flow. It follows the design steps in the traditional design flow, 
but feeds the layout result in the current iteration back to the logic synthesis tools for 
improving the synthesis resuhs in the next iteration. To make such a "construct-by- 
correction" approach effective, the correction step need to be done incrementally with 
respect to the information fed back by layout. However, a different approach 
completely remaps the entire circuit based on the information fed back from the layout 
design, making it difficult to guarantee any convergence when performing the iteration 
between layout and synthesis. 

Another conventional approach is to use a concurrent design flow, which 
performs logic synthesis/technology mapping and placement/routing concurrently. 
However, the optimality of such an approach usually holds for very special circuit 
structures (such as trees) and the main difficulty associated with this approach is its high 
computational complexity. 

Clearly, a better technique is needed for an effective and efficient layout-driven 
synthesis flow. Such a technique should consider layout information during synthesis 
and design planning, such as by combining logic partitioning with retiming and proper 
consideration of global and local interconnect delays, or by exploiting fast interconnects 
available in many PLD architectures during technology mapping. 

As the capacity of PLD devices increases, hierarchical architectures are being 
more widely used, where basic programmable logic blocks, such as look-up tables 
(LUTs) or macrocells, are grouped into a logic cluster and connected by local 
programmable interconnects inside the cluster. There are basically two types of 
clusters, hard-wired connection-based clusters (HCC) and programmable interconnect- 
based clusters (PIC). The layout-driven synthesis flow described in the Layout-Driven 



Timing Optimization section below is mainly targeted for the PIC-based FPGA 
architectures, although, in other embodiments, other architectures are targeted. 

In a PIC, a group of basic logic blocks are connected by a local programmable 
interconnection array that usually provides full connectivity and is much faster than 
5 global or semi-global programmable interconnects, A number of commercial PLDs use 

the PIC architecture, such as the logic array block (LAB) in Altera FLEX lOK and 
APEX 20K devices, and the MegaLAB in APEX 20K devices. For example, in FLEX 
lOK devices (see Figure 1), each LAB consists of eight 4-LUTs connected by the local 
interconnect array. Multi-level hierarchy can be formed easily using PICs, in which a 

10 group of small (lower-level) PICs may be cormected through a programmable 
interconnect array at this level to form a larger (higher-level) PIC. For example, in 
Altera APEX 20K FPGAs (see Figure 2), each LAB consists of ten 4-LUTs connected 
by local interconnects, which forms the first-level PIC. Then, 16 such LABs, together 
with one embedded system block and another level of programmable interconnects, 

15 form a second level PIC, called MegaLAB. Finally, global interconnects are used to 

route between MegaLAB structures and to I/O pins. 

Conventional PLD synthesis algorithms often transform a given design into a 
flat netlist of basic programmable logic blocks (such as LUTs or macrocells) without 
consideration of the device hierarchy. Therefore, a substantial challenge in this area is 

20 to be able to synthesize a given design directly into a multi-level hierarchical 

architecture, with consideration of different interconnect delays and clustering 
constraints at each level. 

Summarv of the Invention 

25 Aspects of the present invention provide a solution for the timing-driven logic 

resynthesis problem, which provides for faster logic synthesis with improved resuUs as 
compared to conventional synthesis techniques. Further, the disclosed timing-driven 
logic synthesis technique may be applied during a technology independent stage after 
area oriented logic optimization is performed to minimize or reduce the depth of the 

30 circuit, or after technology mapping is performed. The timing-driven logic resynthesis 
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techniques can also be integrated into a layout-driven synthesis flow to reduce the 
overall circuit delay. 

One embodiment of the novel timing-driven logic resynthesis technique includes 
(i) a methodology for the general timing-driven iterative refinement-based approach, 
5 (ii) a novel method, named TDO (timing-driven optimization), for optimizing the circuit 

depth after area oriented logic optimization is performed, and (iii) a layout-driven 
synthesis flow that integrates performance-driven technology mapping and clustering 
with TDO to account for the effect of mapping and clustering during the timing 
optimization procedure of TDO. 
10 Embodiments incorporating the novel methodology for the general timing- 

driven iterative refinement-based approach have some or all of the following 
characteristics and advantages: 

1. The need for the (incremental) timing analysis during the iterative 
refinement procedure is reduced or eliminated. 
15 2. The local transformation is able to see or determine more accurate timing 

information (e.g., the arrival times at the transitive fanin signals) so that 
the transformations can be conducted in a more meaningful way to 
reduce the circuit delay. 
3. Design preferences can be much more easily considered (e.g., in hybrid 
20 FPGAs with both LUT clusters and Pterm blocks, it is better to use the 

same logic resources consecutively on a critical path so that the 
subsequent clustering procedure can pack these implementations into one 
cluster to reduce the circuit delay) because of the flexibility of the 
methodology. 

25 4. A general framework is provided which allows the integration of several 

types of local transformations, such as logic resynthesis, mapping, 
clustering, and so on, to enable an integration of currently separated 
design processes. 

The TDO method integrates the novel methodology for the general timing- 
30 driven iterative refinement-based approach and the area recovery technique using 
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restrictive iterative resubstitution. It is able to outperform the state-of-the-art algorithms 
consistently while significantly reducing the run time. 

In one aspect of the present invention, there is a . . . 
[Ray will complete this section based on the independent claims] 

Brief Description of the Drawings 
Figure 1 is a device block diagram of an exemplary programmable logic device 

(PLD). 

Figure 2 is a block diagram of a logic array block structure of another PLD. 

Figure 3 is a diagram of where the timing-driven logic resynthesis aspect of the 
invention may be applied in an exemplary design flow. 

Figure 4 is a flow chart of a generic iterative refinement timing optimization 
procedure {oldTimingOptimize), 

Figure 5 is a diagram showing an exemplary local transformation where a local 
region is resynthesized to reduce logic delay. 

Figure 6 is a diagram showing a deficiency of the oldTimingOptimize procedure 
where an exemplary local transformation does not use accurate timing information. 

Figure 7 is a flow chart of a novel iterative refinement timing optimization 
process {newTimingOptimize). 

Figure 8 is a flow chart of the recursive delay reduction process (reduceDelay) 
shown in Figure 7. 

Figure 9 is a diagram showing exemplary results of executing the recursive 
delay reduction process reduceDelay shown in Figure 8. 

Figure 10 is a flow chart of the local transformation process (transform) shown 
in Figure 8. 

Figure 1 1 is a graph showing the impact of timing optimization on clustering 
and final layout design in relation to area-oriented logic optimization. 

Figure 12 is a flow chart of a novel layout-driven timing optimization process 
{layoutDriven TDO). 




Detailed Description of the Preferred Embodiments 
The following detailed description presents a description of certain specific 
embodiments of the present invention. However, the present invention may be embodied 
in a multitude of different ways as defined and covered by the claims. In this description, 
5 reference is made to the drawings wherein like parts are designated with like numerals 

throughout. 

A novel methodology or process of timing optimization based on iterative 
refinement will be described. This novel methodology reduces or eliminates the need 
for incremental timing analysis during the iterative refinement procedure and any local 

10 transformation is able to utilize more accurate timing information. The timing 

optimization methodology utilizes a delay reduction process to reduce the delay of node 
V. The delay reduction process attempts to recursively reduce the delay of the critical 
fanins of node v instead of conducting the local transformation for v directly. 
Furthermore, in one embodiment, the fanins of node v are sorted in non-ascending order 

15 according to their slack values. Thus, the fanins that have bigger negative slack values 

and are easier to speed up are processed before those fanins that have smaller negative 
slack values and are more difficult to speed up. The novel optimization methodology is 
able to outperform the .state-of-the-art algorithms consistently while significantly 
reducing the run time. 

20 In an exemplary design flow 300 shown in Figure 3, a timing-driven logic 

synthesis technique 310 may be applied during a technology independent stage after the 
area oriented logic optimization 320 is performed to minimize the depth of the circuit, 
or after technology mapping 330 is performed. The timing-driven logic synthesis 
technique 310 can also be integrated into a layout-driven synthesis flow (e.g., after 

25 circuit clustering 340, after placement 350, or after routing 360) to reduce the overall 

circuit delay. 

The remainder of this document is organized as follows: a Problem 
Formulations and Preliminaries section, a Timing-Driven Logic Optimization section, a 
Layout-Driven Timing Optimization section, an Experimental Results and Comparative 
30 Studies section, and a Conclusions section. The Timing-Driven Logic Optimization 

section discusses a novel methodology for performance-driven iterative refinement- 
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based approaches and a novel timing-driven optimization method. The impact of 
technology independent timing optimization on circuit performance after technology 
mapping, clustering and layout design is also discussed. 



PROBLEM FORMULATION AND CONCEPTS 
A Boolean netv^ork N may be represented as a directed acyclic graph (DAG) 
where each node represents a logic gate, and a directed edge <ij> exists if the output of 
gate / is an input of gate j. Primary input (PI) nodes have no incoming edge and 
primary output (PO) nodes have no outgoing edge. Input(v) is used to denote the set of 
fanins of gate v, and output(v) is used to denote the set of nodes which are fanouts of 
gate V. Given a subgraph H of the Boolean network, input(H) denotes the set of distinct 
nodes outside H which supply inputs to the gates in //. Node u is the transitive fanin or 
predecessor of node v if there is a path from u to v. Similarly, node u is the transitive 
fanout of node v if there is a path from v to w. The level of a node v is the length of the 
longest path from any PI node to v. The level of a PI node is zero. The depth of a 
network is the highest node level in the network. A Boolean network is K-bounded if 
\input(v)\ < K for each node v in the network. 

For a node v in the network, a fanin cone (also referred to as a predecessor cone 
or transitive fanin cone) at v, denoted Cy, is a subgraph consisting of v and its 
predecessors such that any path connecting a node in Cy and v lies entirely in Cy. The 
root of Cv is called v. Cy is K-feasible if \input(Cy)\ < K. For. a node v in the network, a 
fanout cone (also referred to as a transitive fanout cone) at v, denoted Dy, is a subgraph 
consisting of v and its transitive fanouts such that any path connecting v and a node in 
Dy lies entirely in Dy. The root of Dy is called v. 

The delay modeling of digital circuits is a complex issue. Certain basic 
concepts, well known to one of ordinary skill in the art, are presented herein to more 
fixlly illustrate the operation of the algorithms. In a Boolean network, it is assumed that 
each node has a single output with possibly multiple fanouts and zero or one edge or 
interconnect from one node to another node. The timing concept for single-output 
nodes can be easily generalized to the case of multiple-output nodes. The concepts of 
pin-to-pin delay and edge delay are defined as follows. 
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Definition 1: The pin-to-pin delay di(v) of a node v in a Boolean Network N is 
the propagation delay from the ith input (pin) of node v to the output (pin) of node v. 

Definition 2: The edge delay d(u,v) of a edge <w,v> in a Boolean network N is 
the propagation delay from the output (pin) of the node u to the corresponding input 
5 (pin) of the node v. 

In the delay models where the interconnection (edge) delay is a constant, the 
edge delay may be combined into the pin-to-pin delay for a more simplistic, yet 
sufficiently accurate, delay modeling. 

Given the propagation delays of each node and connections in a netlist, each PI, 
10 PO or output of every node v is associated with a value called the arrival time t(v), at 
which the signal it generates would settle. The arrival times of the primary inputs 
denote when they are stable, and so the arrival times represent the reference points for 
the delay computation in the circuit. Often the arrival times of the primary inputs are 
zero. Nevertheless, positive input arrival times may be usefiil to model a variety of 
15 effects in a circuit, including specific delays through the input pads or circuit blocks that 

are not part of the current logic network abstraction. 

The arrival time computation may be performed in a variety of ways. A model 
is considered herein that optionally divorces the circuit topology from the logic domain, 
i,e,, arrival times are computed by considering the dependencies of the logic network 
20 graph only and excluding the possibility that some paths would never propagate events 
due to the specific local Booleaii functions. The arrival time t(v) at the output of each 
node V may be computed as follows. Let w/ be the zth fanin node of v, then 

t(v) = max (t(ui) + d(uiv) + di(v)) (1) 

25 0<i<\input(v)\ 

The arrival times may be computed by a forward traversal of the logic network 
in 0(n+m) time, where n and m are the number of nodes and number of edges in the 
network, respectively. The maximum arrival time occurs at a primary output, and it is 
30 called the critical delay of the network. The propagation paths that cause the critical 

delay are called critical paths. 
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The required time at the output of every node v, denoted as i (v), is the required 
arrival time at v in order to meet the overall circuit timing constraint as deflned by 

{Songjie, please fill in}. The required times may 

be propagated backwards, from the POs to the Pis, by way of a backward network 
traversal. Let be the zth fanout node of v, and v be the yth fanin of node w^, then: 

t (v) - min (i (uj) - dj(u0 - d(v, uj)) (2) 

o</<|outpw/('v;i 

The difference between the required time and the actual arrival time at the output 
of each node is referred to as timing slack, namely: 

s{v)=t(v)-t(v) (3) 

The required times and the timing slacks may be computed by the backward 
network traversal in O(n^m) time, where n and m are the number of nodes and number 
of edges in the network, respectively. Critical paths are identified by nodes with zero 
slack, when the required times at the primary outputs are set equal to the maximum 
arrival time. 

For each node v in the Boolean network A^, d-fanin-cone and d-critical-fanin- 
cone may be defined as follows. 

Definition 3: The d-fanin-cone of a node v in a Boolean network N is the set of 
nodes that meet the following requirements: 

1 . They are in the transitive fanin cone of node v; and 

2. they are at most distance d (d levels of logic) away from the node v. 
Definition 4: The d-critical-fanin-cone of a node v in a Boolean network N is 

the set of nodes that meet the following requirements: 

1 . They are in the d-fanin-cone of node v; and 

2. each of them has a negative slack. 

The timing optimization problem for multi-level networks during the technology 
independent stage may be formulated as follows: 
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Problem 1: Given a AT-bounded Boolean network N, transform 
to an equivalent /T-bounded network so that the circuit depth is 
minimized. 

The general timing optimization problem for multi-level networks concerning 
5 the circuit delay after the layout design may be formulated as follows: 

Problem 2: Given a AT-bounded Boolean network N, transform 
to an equivalent A^-bounded network A^' so that the circuit delay after the 
layout design is minimized. 

10 TIMING-DRIVEN LOGIC OPTIMIZATION 

Conventional timing-driven logic optimization is typically conducted using an 
iterative refmement-based approach on the critical netlist due to the potential area 
penalty associated with the local transformations for delay reduction. In each pass (or 
iteration), the overall circuit delay target, which is smaller than the current maximum 

15 arrival time is set, and a set of local transformations are selected to meet this delay 

target. An exemplary local transformation is shown in Figure 5 where a local region is 
resynthesized to reduce delay. Node v (510) before resynthesis has t=4, while after 
resynthesis, node v (510') has t=2, for a delay savings of 2. If the delay target is met, 
the whole process is disadvantageously repeated to meet a new delay target. The timing 

20 optimization stops when it is not possible to reduce the circuit delay any more. In 
contrast to the conventional method, one embodiment of the present invention 
advantageously reuses delay optimizations from each iteration. 

The novel methodology for the general timing-driven iterative refmement-based 
approach is described in subsection A below. This methodology is applied to the novel 

25 method, called herein TDO (timing-driven optimization), for optimizing the circuit 

depth after the area oriented logic optimization in subsection B. The impact of 
technology independent timing optimization on the circuit performance after technology 
mapping, clustering and layout design is discussed in subsection C. 



-13- 



A. A Novel Methodology for Performance-Driven Iterative Refinement Approaches 
The present invention, which utilizes performance-driven iterative refinement 
based on timing considerations, offers several important advantages over conventional 
approaches. Conventional timing optimization algorithms generally use the generic 
iterative refinement procedure 400 shown in Figure 4. One embodiment of the 
procedure 400 may be performed by the pseudo-code shovra in Table 1. 

This procedure template 400 {oldTimingOptimize) may be customized to yield 
specific algorithms by changing state 404 and state 406 of the optimization loop, or by 
using different transformations at state 408. The oldTimingOptimize procedure 400 was 
discussed in (Singh, KJ., Performance Optimization of Digital Circuits, Ph.D, 
Dissertation, University of California at Berkeley^ 1992), where each of the three states 
404, 406, 408 was studied and the appropriate strategies applied. As the result, the 
algorithm proposed in (Singh, K.J., Performance Optimization of Digital Circuits, Ph.D. 
Dissertation, University of California at Berkeley, 1992) is able to generate solutions 
with better qualities than the previous approaches using iterative refinement procedures, 
however, with a rather long and unpredictable computation time. Both the good quality 
of the solution and the long run time are due to the Binary Decision Diagram (BDD) 
based approach used in state 406 (transformation selection). A recent study (Tamiya, 
Y., Performance Optimization Using Separator Sets, in Proa. Intl Conf On Computer 
Aided Design, pages 191-194, 1999) tries to speed up the transformation selection 
procedure by computing multiple separator sets instead of using BDDs. 



TABLET 



Template oldrimingOptimize(N) 



3 



2 



1 



repeat 

set the delay target and determine the critical region 
select a set of critical nodes to be transformed 
apply the transformations on the selected nodes 



until the delay cannot be reduced or constraints are violated 



traditional timing optimization procedure 
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Instead of perfecting the oldTimingOptimize flow, it was discovered, however, 
that this procedure has the following intrinsic disadvantages: 

1 . When the local transformations are evaluated for delay reductions at each 
5 critical node at state 406, it is assumed implicitly that the arrival times of 

the transitive fanins of this node will not be changed. This is generally 
not true with the iterative refinement-based approach. Thus, the local 
transformations conducted on the selected nodes at state 408 do not use 
the most up-to-date timing information, such as seen in Figure 6. While 

10 the algorithm in (Pan, P., Performance-Driven Integration of Retiming 

and Resynthesis, in Proc. Design Automation Conf.^ pages 243-246, 
1999) does allow the transformations to use the accurate timing 
information by propagating the arrival times fi-om Pis to POs, it 
disadvantageous^ invalidates the assumptions for the node selection at 

1 5 state 406, which are based on the old timing information. 

2. Depending on how the transformation at state 408 is conducted (whether 
to propagate the arrival times from Pis to POs), it may require the 
(incremental) timing analysis after each iteration of the delay reduction. 

Overcoming these problems is crucial for the timing optimization of today's 
20 high-performance FPGAs. Therefore, there is a need for new methodology in 
performing timing optimization based on the iterative refinement. One embodiment of 
the novel methodology is described in conjunction with Figure 7 and Figure 8. 

The novel methodology or process of timing optimization based on iterative 
refinement eliminates or reduces the need for incremental timing analysis during the 
25 iterative refinement procedure. The novel methodology enables any local 

transformation to be able to utilize more accurate arrival times at the transitive fanin 
signals. The timing optimization methodology utilizes a delay reduction process to 
reduce the delay of node v. The delay reduction process attempts to recursively reduce 
the delay of the critical fanins of node v instead of conducting the local transformation 
30 for V directly. The critical PO nodes may sorted according to their slack. The PO nodes 

that have bigger negative slack values and, thus, are easier to speed up, are processed 
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before the PO nodes that have smaller negative slack values, as will be further described 
below. 

Referring to Figure 7, the overall flow of a newTimingOptimize process 700 will 
be described. Portions of states 704 and 706 (described hereinbelow) may be 
customized by the specific timing optimization method. One embodiment of the 
process 700 may be performed by the pseudo-code shown in Table II. 



TABLE II 





procedure fi€.vj7%Tfii7iQ0ptiTTi%Z€>{^N') 


1 


do an initial timing analysis with arrival times computed for every node in N" 




and circuit^elay is the maximum arrival time in the circuit 


2 


set the initial delay reduction step reduce^tep* 


3 


repeat 


4 


cki^elayUarget = ckiulelayr'reducestep 


5 


success — true ■ 


6 


for eadi critical primary output .^^ in non- ascending or 




according to the current negative slack do 


7 


i( reducebelay(Vi ckiuielaydarget) —— / then 


8 


success — false ■ • >^ .' ' ■■■ = ; . 


9 


■ break • • 


10 


\f success true th^n . 


11 


cktuieldy = cktudelayJarget 


12 


else 


13 


njd'}nRt reduces tep* 


14 


prepare for the next iteration* 


15 


until the delay cannot be reduced or constraints are violated 



Beginning at a start state 702, process 700 proceeds to state 704 and performs an 
initial timing analysis. The initial timing analysis computes the arrival times for each 
node in N. The circuit delay, denoted as ckt_delay, is the maximum arrival time. An 
initial delay reduction step is then selected, which controls the pace of the timing 
optimization. Continuing at state 706, the first operation in each timing optimization 
iteration of process 700 is to choose or set the delay target {ckt delay Jar get) based on 
the current circuit delay (cktjdelay) and a given delay reduction step {reduce step). A 
flag labeled success indicates whether the delay target may be met in one iteration of the 
timing optimization and it is set initially to true. Based on the delay target, each critical 
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PO has a negative slack. In one embodiment, these critical POs are sorted in non- 
ascending order according to their slack. This order ensures that the POs that have 
bigger negative slack values and, thus, are easier to speed up, are processed before those 
POs that have smaller negative slack values and, thus, are more difficult to speed up. 
5 For example, if there are three critical POs with slack values of -2, -1 and -3, 

respectively, then the PO with the slack of -1 is processed first, and the PO with the 
slack of -3 is the last to be processed. The rationale for this ordering is that after 
processing the POs that are easier to speed up by local transformations, the delay 
savings resulting fi"om those transformations may be used or shared by the delay 

10 optimization of those POs that are more difficult to speed up so that the POs with 
smaller negative slacks may be sped up in the same iteration as those POs with bigger 
negative slack values. This strategy is very effective and permits a big delay reduction 
step {reduce step) to be used, which results in both area savings and delay reductions. 

Proceeding to process 708, for each critical PO v and its delay target, a recursive 

15 delay reduction (reduceDelay) procedure is invoked. Process 708 will be fiirther 

described in conjunction with Figure 8. If the delay of the critical PO V is successfully 
reduced by local transformations in its transitive fanins, as determined at a decision state 
710, process 700 proceeds back to state 706 to operate on the next critical PO. 
Otherwise, the flag success is set to false and there is no need to continue the delay 

20 reduction for the remaining POs. 

Continuing at state 706, if the delay targets of all the critical POs are met, the 
circuit delay (ckt_delay) is set to the current delay target ckt delay target and the next 
iteration of timing optimization may begin. If the delay target is not met for one or 
more POs, the delay reduction step (reducejstep) may be adjusted to a less aggressive 

25 value, which may be done in a customized manner. {Songjie, please explain how it 

may be customized} Preparation for the next iteration may then be performed, possibly 
including a recovery to the previous netlist without conducting the partial 
transformations. 

Referring to Figure 8, the recursive delay reduction (reduceDelay) process 708 
30 for a node v with respect to a specific delay target delay target will now be described. 

Portions of process 810 (described hereinbelow) may be customized by a specific 
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timing optimization method, such as timing-driven decomposition, timing-driven 
cofactoring, generalized bypass transform, or timing-driven simpHfication. One 
embodiment of the process 708 may be performed by the pseudo-code shown in 
Table III. 



TABLE m 





procedure reduceDelay{v, delay Jarijet) 


1 


update the arrival time t(v) of v according to its fanins' arrival times 




if t{v) < delay.target then 


3 


return true 


4 


if t; is a primary input then 


5 


return false 


6 


for each fanin n € input(v) in non-ascending order 




according to the current slack do 


7 


fanin-delayJtarget — delay -tar get - di(v) — rl(ii, v) 




y* n is v's ith fanin node */ 


8 


if t{u) > fanin^elayJtar get then 


, 9 


if r educe Delay{Ui fanirudelayJ^argei) —— ■false then 


10 


return transform{Vi delay Jar get)* 




/* transform returns true if delay-target is met, otherwise, 




it returns false */ 


11 


update the arrival time t(v) of v according to its fanins' arrival times 


12 


return true 



recursive delay reduction for a node with respect to a delay target 



Beginning at a start state 802, process 708 moves to state 804 and updates the 
arrival time t(v) of node v according to its fanins' arrival times, corresponding edge 
delays and pin-to-pin delays (refer to Problem Formulation section above). This update 
is performed as some of the fanins of node v may have been sped up during the timing 
optimization on other critical paths. If the updated arrival time already meets the delay 
target, then a Boolean true is returned. If v is a primary input, then a Boolean false is 
returned indicating that v cannot speed up. Otherwise, the fanins of node v are sorted in 
non-ascending order according to their slack values, in one embodiment. This order 
ensures that the fanins that have bigger negative slack values and are easier to speed up 
are processed before those fanins that have smaller negative slack values and are more 
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difficult to speed up. For example, if v has two fanins uj and U2 with slack of -2, -1, 
respectively. Then the fanin u j with slack of -1 is processed first, and the fanin U2 with 
slack of -3 is the last to be processed. The rationale for this ordering is similar to that of 
the critical PO ordering used in the newTimingOptimze process 700. 

A feature of the reduceDelay process 708 is that in order to reduce the delay for 
node V, instead of conducting the local transformation for v directly, the process 
recursively reduces the delay of the critical fanins of v, where possible. It is clear 
{Songjie, please explain why} that in this procedure, the transformations closer to the 
primary inputs are preferred. 

For each fanin u of node v, its delay target (fanin jielayjqrget) is computed 
according to the delay target for v, the pin-to-pin delay di (assuming that u is the rth 
fanin of node v) and the edge delay d(u,v). For a fanin u with t(u) > fanin delay Jar get, 
a reduceDelay process 708' is invoked on u for a recursive delay reduction. Continuing 
at a decision state 808, if for some fanin u of v, the delay reduction is not successful, 
then the reduceDelay process 708 stops trying to reduce the delay for other critical 
fanins of v. Instead, process 708 proceeds to a transform process 810 to conduct the 
local transformation on the node v itself, with the goal of hitting the delay target. 

An example of recursive delay reduction is shown in Figure 9. The arrival times 
t=7 at fan-in u3, t=8 at ul, t=8 at u6, t=9 at u2 and t=10 at v represent initial arrival 
times before delay reduction for the exemplary logic circuit. After delay reduction, the 
final arrival times are t=5 at u3, t=5 at u4, t=6 at ul, t=7 at u2 and t=8 at v. After 
applying a local transformation at node u2 as part of the delay reduction process, u2 no 
longer depends on u6 (u6 is not the fanin of u2 any more) and u2's arrival time is 
reduced from 9 to 7 so that the arrival time of v is 7+1=8. 

If it turns out that all the critical fanins of node v can be sped up to meet their 
delay targets as determined at the decision state 808, then there is no need to conduct 
any transformation on v itself Instead, the arrival time t(v) of v is directly updated and a 
Boolean true is returned indicating that the delay reduction for node v with respect to 
the delay target was successfiil. 

From the above description of the novel methodology of the timing optimization 
based on the iterative refinement, it will be understood by one of ordinary skill in the art 
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that it is general enough to consider different pin-to-pin delays and distinctive edge 
delays. This general procedure may be customized for any time-driven optimization 
that adopts the iterative refinement-based approach. In the next subsection, this 
procedure is applied to the timing optimization during the technology independent 
stage. The experimental resuhs show that this new method produces very favorable 
results compared to conventional algorithms. 

The novel methodology for the general timing-driven iterative refinement-based 
approach has the following advantages, though particular embodiments may include 
only some of the advantages: 

1 . It reduces or completely eliminates the need for the (incremental) timing 
analysis during the iterative refinement procedure. 

2. It allows the local transformation to be able to see more accurate timing 
information (e.g., the arrival times at the transitive fanin signals) so that 
the transformations may be conducted in a more meaningful way to 
reduce the circuit delay. 

3. Its flexibility makes it much easier to consider the design preferences 
(e.g., in hybrid FPGAs with both LUT clusters and Pterm blocks, it is 
better to use the same logic resources consecutively on a critical path so 
that the subsequent clustering procedure may pack these implementations 
into one cluster to reduce the circuit delay.) 

4. It provides a general fi-amework to integrate several types of local 
transformations, such as logic resynthesis, mapping, clustering, and so 
on, to enable an integration of currently separate design processes. 

B. A Novel Method for Timing-Driven Optimization 

Based on the timing optimization fi-amework presented in the previous 
subsection, the novel method, termed herein TDO (timing-driven optimization), for 
optimizing the circuit depth after the area oriented logic optimization will now be 
fiirther described. In one embodiment, the input to TDO is a 2-bounded netlist and the 
output of TDO is also a 2-bound netlist. The novel method may be obtained by (i) 
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# # 

customizing portions of states 704 and 706 (Figure 7) and process 810 (Figure 8), and 
(ii) performing the area recovery after each successfiil circuit delay reduction. 

The customization involves portions of states 704 and 706 of the 
newTimingOptimize process 700 (Figure 7) and process 810 of the reduceDelay process 
5 708 (Figure 8). The customization of state 704 in newTimingOptimize process 700 is to 

determine the initial delay reduction step reduce step, which is adjusted to a less 
aggressive value if a failure occurs in meeting the delay target based on the current 
reduction step. The value reduce _step may be any number fi-om one to the overall 
desired delay reduction that is the difference between the initial circuit delay and the 
10 ultimate delay target. If reduce _step is too small, for example, one (1), the timing 
optimization may proceed in a very slow fashion, which has several drawbacks: 

• The critical path information during one iteration of the circuit delay 
reduction is rather limited as it does not consider the forthcoming critical 
paths after the current iteration, which, if explored together with the 

15 current critical paths, may yield more optimal results for the area and 

also the delay in the long run of the timing optimization. 

• The overall timing optimization time would be long. 

However, if reduce_step is too large, the timing optimization works on a large 
set of critical nodes. However, it is unlikely that the delay target can be achieved, in 

20 which case reduce_step has to be adjusted to a less aggressive value and the whole 
procedure has to be restarted. One can use the well-known method described in (Singh, 
K.J., Performance Optimization of Digital Circuits, Ph.D. Dissertation, University of 
California at Berkeley, 1992) to compute the lower bound of the delay reduction at the 
beginning of every iteration and use that as reduce_step. However, computing the 

25 lower bound of the delay reduction involves conducting the transformations on every 

critical node, which may be a timing consuming process, especially when the timing 
optimization approaches the end, where more nodes become critical. In the 
experimental results, an empirical value of four (logic levels) is chosen beforehand as 
the initial reduce step, 

30 The customization of state 706 in the newTimingOptimize process 700 involves 

adjusting reduce step to a less aggressive value if another failure occurs in meeting the 
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delay target based on the current reduction step. This may be accomplished by 
decrementing reduce_step by one or some other predetermined value. The resultant 
reduce jstep may be used as the delay reduction step until a failure occurs in meeting the 
delay target. 

5 Further customization of state 706 in the newTimingOptimize process 700 is 

performed if the delay target cannot be achieved for some PO. If so, the previous netlist 
is advantageously recovered as the starting point to the next delay reduction iteration 
based on a less aggressive reduce step without conducting the partial transformations 
(transformations may been conducted on some nodes to reduce their delays.). 
10 The customization of process 810 in the reduceDelay process 708 involves 

determining the transformation method, Le,, which type of transformation is to be 
performed to reduce the delay of node v to meet its delay target. The transformations 
that alter the structure of a part of the circuit, such that the delay through the part is 
reduced, include, but are not limited to: 
15 1. Timing-driven decomposition. As is well known in the art, timing- 

driven decomposition decomposes a complex function / into a Boolean 
network composed of 2-input functions having minimum delays. This is 
done primarily through the extraction of divisors that are good for 
timing. Whether a divisor is good or not depends on the arrival times of 
20 its supporting inputs and the area saving it may provide. A good divisor 

should not have late arriving signals as its inputs. The best divisor g is 
chosen each time and substituted into /. The function g is then 
decomposed recursively, followed by the recursive decomposition of the 
resulting function /. The choice of divisors that are evaluated affects the 
25 quality of the decomposition. Algebraic divisors, such as kernels or 

two-cube divisors, are well-known techniques that may be used for the 
extraction. If, after a predetermined number of attempts, the dividend 
function v does not have any good divisors, v is a sum of disjoint-support 
product terms and its decomposition into 2-input functions may be 
30 performed using a conventional Huffinan-tree based structural 

decomposition procedure. 
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2. Timing-driven cofactoring. The timing-driven cofactoring technique is 
a well-known technique for performing optimization. Given a function 
/, the latest arriving input x is determined. / is then decomposed as / = 
^fx^ ^'fx* (6) Ux is / with the input x set to one, and /x' is / with x set 
5 to zero). A straightforward implementation realizes fx and fx* 

independently, which may result in a large area overhead. The overhead 
may be reduced if logic sharing between fx and fx* is considered. This 
technique is a generalization of the design of a carry-select adder 
technique. 

10 3. Generalized bypass transform. The basic idea in this method is to 

change the structure of the circuit in such a way that transitions do not 
propagate along the long paths in the circuit. Given a function / and a 
late arriving input x, g ^ fx® fx' represents conditions under which / 
depends on x. g is then used as the select signal of a multiplexer, whose 

15 output is /. If g is one, the output / is simply x or x\ If g is zero, the 

output is the function / with x set to either zero or one. The long path 
that depended on x is replaced by the slower of the two functions: g and 
g with X set to a constant. This transformation is a generalization of the 
technique used in a carry-bypass adder. 

20 4. Timing-driven simpliflcation. As discussed above, this simplification 

computes a smaller representation of a function using a don't care set 
that is derived from the network structure and also possibly from the 
external environment. The goal of timing-driven simplification is to 
compute a representation that leads to a smaller delay implementation. 

25 This may be achieved by removing late arriving signals from the current 

representation using appropriate don't care minterms, and substituting 
therefore early arriving signals. 
The transformation based on the timing-driven decomposition may generally 
produce the best resuhs in terms of the delay reduction among the transformations listed 

30 above. Therefore, the timing-driven decomposition is used in process 810 method with 

kernels as possible divisors. 
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Based on the above discussion, the customization of process 810 (transform(v, 
delay Jar get) in reduceDelay (Figure 8) is shown in Figure 10. One embodiment of the 
process 810 may be performed by the pseudo-code shown in Table IV. 



TABLE IV 





procedure trans fofTn{v^ delay JLorget) 


1 


update the arrivsJ time t{v) of v according to its fanins' arrival times 


2 


\f t(ii\ <C d^loAt -inr Q€.t then 


3 


return tvue 


4 


if v is a primary input then 


5 


return false 


6 


set the minimum collapse depth min-d and maximum collapse depth maxjd 


7 


d = minjd 


8 


while d < maxjd do 


9 


collapseCTitial{v^d) 


10 


timing Decompose(v) 


11 


if t(v) < delay Jar get then 


12 


return true 


13 


collapse{vtd) 


14 


timing Decompose{v) 


15 


if t{v) < delay -tar get then 


16 


return true 


17 


d-d-\-\ 


18 


return false 



delay reduction for a node with respect to a delay target 



Beginning at start state 1002, process 810 moves to state 1004 and updates the 
arrival time of node v. Proceeding to state 1006, in order to apply the timing-driven 
decomposition on node v, a transformation region that is a partial fanin cone of v is 
formed and collapsed into v. As in most of the other timing optimization algorithms, 
the collapse depth d is used to control the size of the transformation region. The choice 
of the collapse depth d certainly influences the quality of the TDO method. A large d is 
useful in making relatively large changes in the delay since a larger region that results in 
a more complex function provides greater flexibility in restructuring the logic. 
However, this results in both longer run time, due to the collapse operation and the 
timing-driven decomposition procedure, and a bigger area overhead as more logic is 
duplicated. An empirical value of three has been chosen for d in the past. 
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Experimentation results show that when d is larger than three, the area overhead may be 
unwieldy. Therefore, in one embodiment, the maximum collapse depth is set to three, 
though other values may be used as well. If a smaller collapse depth may help meet the 
delay target, it can be used to reduce the area overhead. In one embodiment of the TDO 
5 method implementation, a value of two is used as the minimum collapse depth. 

Advancing to state 1008, a variable d is set to the value of min-rf, e.g., two. 

Given a certain collapse depth d, either the d-fanin-cone or the d-critical-fanin- 
cone, a subset of d-fanin-cone, may be used. Using the d-fanin-cone generally results in 
a better delay reduction as compared to using the d-critical-fanin-cone, with, however, a 
10 larger area overhead. In one embodiment, as the overall TDO method is run time 

efficient, the method tries to reduce delay using the d-critical-fanin-cone first, and if 
that fails, the d-fanin-cone is used. Proceeding to a decision state 1010, process 810 
Q determines if d is less than or equal to the value of mzx-d, e.g., three. If so, process 810 

SI moves to collapseCritical process 1012. The collapseCritical process 1012 collapses 

J 15 the d-critical-fanin-cone for v based on depth d. Collapsing a sub-netlist N' means to 

£ eliminate all the internal nodes of N' (but not the output nodes of N'). This is a well- 

ffl known logic operation in multi-level logic optimization. 

□ Continuing at a timingDecompose process 1014, the timing-driven 

^ decomposition is performed on v. Advancing to a decision state 1016 process 810 

W 20 determines if the delay target is met. If so, the local transformation process 810 
M: completes and returns with a true condition at state 1028. If the delay target is not met, 

as determined at decision state 1016, process 810 moves to a collapse process 1018. 
Process 1018 collapses the d-fanin-cone for v based on depth rf. The collapse process 
1018 is similar to that of the collapseCritical process 1012 except for the fan-in cone 
25 used in the processes. At the completion of process 1018, execution continues at 

timingDecompose process 1014', which is similar to process 1014 described above. If 
the delay target is met, as determined at a decision state 1022, process 810 completes 
and returns with a true condition at state 1028. If the delay target is not met, process 
810 proceeds to state 1024 and increments depth d by one and moves back to decision 
30 state 1010 as described above. If d is determined to be greater than max- J at decision 
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state 1010, process 810 completes without meeting the delay target and returns with a 
false condition at state 1026. 

Another feature of TDO includes performing an effective area recovery after 
each successful circuit delay reduction. Each delay reduction iteration (states 706 to 
710 (Figure 7)) involves the transformations on a set of critical nodes, which result in 
the duplication of logic due to the collapsing. The area recovery feature removes the 
redundant nodes (2-input nodes for TDO) whose functions or the complements of the 
functions are already present in the network. This may be accomplished by performing 
a restricted resubstitution, i.e., g is resubstituted into / only if f = g or f = g'. A sweep 
operation following the resubstitution may remove the buffers and inverters generated 
by the resubstitution so that the delay at every node in the circuit may not be increased 
while the circuit area is reduced. 

In summary, the TDO method integrates the novel mechanism for the general 
iterative refinement flow and the area recovery technique using the restrictive iterative 
resubstitution. It outperforms the state-of-the-art algorithms consistently while 
significantly reducing the run time. 

C. The Impact of Timing Optimization on Subsequent Design Processes 

With the novel method presented for the timing optimization during the 
technology independent stage, which have proved to be very efficient and generate 
solutions with superior qualities than conventional algorithms, it would be worthwhile 
and interesting to analyze the impact of technology independent timing optimization on 
the circuit performance after technology mapping, clustering and layout design. 

Using the TDO method, the impact of the timing optimization on the subsequent 
technology mapping was analyzed. The operation of Optarea is a technology 
independent area optimization, which is comparable to the area optimization script 
scripLalgebraic in SIS (Sentovich et al., SIS: A System for Sequential Circuit Synthesis, 
Electronics Research Laboratory, Memorandum No. UCB/ERL M92/41, 1992). The 
comparison of the optimization results is based on the resulted 2-bounded netlists. The 
comparison of the mapping results are based on the mapping into a 4-LUT FPGA that is 
performed by the state-of-the-art algorithm that is able to achieve the optimal depth 



-26- 




while minimizing the mapping area (Hwang, Y.-D., Logic Synthesis for Lookup-Table 
Based Field Programmable Gate Arrays, Ph.D, Dissertation; University of California at 
Los Angeles, 1999). The technology independent timing optimization by TDO reduced 
the circuit depth {do) obtained by Optarea by 56% with a 6% area increase (a^^). After 
5 the technology mapping, the delay reduction is decreased from 56% after the 

optimization to 30% {d^n) with an overall 13% increase on area {a^n)^ Of course, in 
other examples, these numbers may vary. 

The impact of the timing optimization on clustering and the final layout design 
was analyzed. The clustering is an optimization step (refer to the Layout-Driven 
10 Synthesis subsection above) to physically group the mapped 4-LUTs into the clusters of 

an FPGA. In this example, it is assumed that each cluster has ten 4-LUTs, which is the 
same as the LAB structure in an APEX 20K device. The comparison of the clustering 
O results are based on the duplication-free clustering performed by the algorithm that is 

SI able to achieve the optimal delay for all the reported benchmarks. The delay after 

J" 15 clustering {dc) is estimated by a timing analysis tool that considers the LUT logic delay, 

£ the intra-cluster interconnection delay and the inter-cluster interconnection delay. The 

3l final layout may be performed by the Quartus version 2000.03 from Altera on the 

p EPF20K400BC652-2 device. Using the present invention, the delay reduction is fiuther 

P decreased from 30% after the mapping to 15% {d^^ after the clustering and 12% after 

fy 20 the final layout design (t//). As long as the mapping is completed, the circuit logic area 
U will not change much, meaning that the packing ratio (the average number of LUTs 

packed into one cluster) achieved by the clustering is more or less a constant. 

Figure 1 1 summarizes the impact of the timing optimization on clustering and 
the final layout design. {Songjie, does this figure show use of the invention?} In 
25 general, if the overall design process is separated into several design optimization 

stages, such as the technology independent optimization, mapping, clustering, and place 
and route, to be performed sequentially, the delay reduction obtained in the earlier 
stages will not be preserved after the optimization by the later stages. The reason has 
been that the optimization done in each stage tends to reduce the delay along the critical 
30 paths resulted from the previous design stages. Therefore, a circuit with a much longer 

critical path depth resulting from the pure area optimization in the technology 
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independent optimization stage gets optimized more than the circuit with a smaller 
depth achieved by the timing optimization in the subsequent mapping, clustering and 
layout design. 

Furthermore, each cluster typically has capacity constraints and also pin 
constraints. Therefore, both circuit depth and area may have an impact on the 
clustering, and ultimately it is the circuit topology that affects the overall clustering 
performance, which is difficult to consider during the technology independent timing 
optimization. 

A graph 1100 shows a comparison between the pure area-oriented logic 
optimization (line 1110) and area-plus-timing optimization (line 1112) on the circuit 
delay after optimization, mapping, clustering and layout design. The Y-axis 1114 
represents circuit delay ratio (discussed below) and X-axis 1116 represents certain 
design stages, including optimization, mapping, clustering and layout. The absolute 
delay values for area-oriented logic optimization are all scaled to one (line 1110), and 
the absolute delay values for area-plus-timing logic optimization are all scaled 
accordingly. As an example of what is meant by delay ratio, 0.44 point on line 1112 
means that with timing optimization, the delay after optimization is only 44% of the 
delay achieved by pure area optimization. 

On the other hand, analysis results suggest that the estimated delay after the 
circuit clustering (dc) correlates reasonably well, in terms of the relativity, with the 
layout delay (df). A clustering-driven synthesis flow is described in the next section, 
which considers the effect of the mapping and clustering during the timing optimization. 

LAYOUT-DRIVEN TIMING OPTIMIZATION 
A layout-driven synthesis flow that considers the effect of technology mapping 
and circuit clustering during the technology independent timing optimization will now 
be described. This layout-driven synthesis flow makes use of mapping and clustering to 
help detect the circuit topology and uses the resulting inter-cluster edges (i.e., the edges 
whose terminals are spread in different clusters) and their delays as guidance for the 
timing optimization procedure. To ensure that the changes made during the timing 
optimization are incremental and will finally converge to a delay reduction after the 
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mapping and clustering, the timing optimization is performed within each cluster. In 
one embodiment, FPGAs with hierarchical interconnection structures, e.g., PIC-based 
FPGAs such as the Altera FLEX lOK and APEX 20K, and the Actel Pro ASIC 500K, 
are targeted. 

Figure 12 shows the overall flow of a layoutDrivenTDO process 1200. One 
embodiment of the process 1200 may be performed by the pseudo-code shown in 
Table V. 



TABLE V 





procedure layoutDrivenTDO{N) 


1 


timingDecompo3e(Ny2) 


2 


oldu^kt-delay = cx> 


3 


repeat 


4 


N' — duplicate(N) 


5 


mapping(N^) 


6 


cluster.ing(N^) 


7 


do a timing analysis on TV' and circuit jdelay is 




the maximum arrival time in the circuit 


8 


if cktAelay < old.cktjielay then 


9 


oldjckiuielay = ckiAelay 


10 


annotate net delays on N based on the clustering results on AT' 


11 


else 


12 


N - TV" is the 2-bounded netlist saved before the last TDO*/ 


13 


suspend the transformations on TV' done in the last TDO 


14 


N" = duplicate(N) 


15 


restrictecH'DOiN) 


16 


until the delay cannot be reduced nr constraints are violated 



layout-driven timing optimization procedure 



Beginning at a start state 1202, process 1200 moves to state 1204 where the 
netlist N is decomposed into a 2-bounded netlist. A variable old_ckt_delay, indicating 
the circuit delay after clustering before the last timing optimization, is initially set to 
infinity. Each timing optimization iteration consists of the operations from state 1206 to 
state 1214. The first operation in each iteration is to dupHcate the 2-boimd netlist N to 
N* at state 1206. A^' may be optimized in the subsequent mapping and clustering 
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procedure. A timing analysis is conducted on the mapped and clustered netlist to 
compute the arrival times for every node in A^'. The circuit delay after clustering, 
denoted as cktjdelay, is the maximum arrival time. 

If the circuit delay (ckt_delay) after the last timing optimization, mapping and 
5 clustering is indeed better than the previous one {old_ckt_delay\ as determined at a 

decision state 1208, the last timing optimization is considered to be good. Thus, 
old_ckt_delay is updated to the current delay and the net delays in the 2-bounded netlist 
N are set up based on the clustering results of A^'. Therefore, the delay model used 
during the timing optimization is the same as the one used by the circuit clustering. 

10 If the circuit delay {cktjdelay) after the last timing optimization, mapping and 

clustering does not improve over the previous one {old ckt delay) as determined at 
decision state 1208, the last timing optimization is considered to be bad. Thus, the 
nethst N is reversed back to the 2-bounded netlist A^" saved before the last timing 
optimization at state 1210 and the transformations done in the last timing optimization 

15 are suspended. Before the next timing optimization, the netUst N is duplicated to A 

restricted timing optimization is then performed on at a restrictedTDO procedure 
1212, where, in one embodiment, the suspended transformations are not considered and 
the optimization is only done within each cluster, which means that any local 
restructuring does not go across the cluster boundary. If the delay cannot be fiirther 

20 reduced or the constraints are violated, as determined at decision state 1214, process 

1 200 completes at an end state 1216. 

At the end of the layoutDrivenTDO(N) process 1200, the netlist A/^ is a 2- 
bounded netlist that has been optimized for delay with the consideration of the potential 
impact on the subsequent mapping and clustering. 

25 In summary, compared with the traditional timing optimization flow, the layout- 

driven synthesis flow that considers the effect of technology mapping and circuit 
clustering during the technology independent timing optimization has the following 
advantages and differences. 

1. The flow layoutDrivenTDO 1200 makes use of the mapping and 

30 clustering to help detect the circuit topology and uses the resulting inter- 

cluster edges and their delays as the guide for the timing optimization 
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procedure. Traditional timing optimization does not consider these 
factors, 

2. In contrast to traditional methods, to ensure that the changes made during 
the timing optimization are incremental and will finally converge to a 
delay reduction after the mapping and clustering, the timing optimization 
in layoutDrivenTDO 1200 is performed within each cluster. 

3. In the procedure restrictedTDO 1212 of the flow, the suspended 
transformations, which may harm or have no evident benefit in reducing 
the circuit delay, are not considered. Traditional timing optimization 
does not suspend any transformations since no information is available 
indicating whether a transformation is good or not. 

4. In the procedure restrictedTDO, the net delays in the 2 -bounded netlist A'^ 
are set up based on the clustering result. Thus, the delay model used 
during the timing optimization is the same as the one used by the circuit 
clustering. The edge delays in traditional timing optimization are zero. 



EXPERIMENTAL RESULTS AND COMPARATIVE STUDY 
To effectively carry out the experimentation, a set of benchmarks are first 
selected. Twenty eight benchmark circuits, which are among the largest in the MCNC 
benchmark suite (Yang, S., Logic Synthesis and Optimization Benchmarks User Guide 
Version 3.0, Technique Report, MCNC, January 1991), are selected for the 
experimentation. 

A comparison was performed between TDO and the SIS speedup algorithm 
(Singh, K.J., Performance Optimization of Digital Circuits, Ph.D. Dissertation, 
University of California at Berkeley, 1992) on the benchmark circuits. On average, the 
solutions generated by speedup have 10% more delay and 14% more area after the 
optimization, and have comparable circuit delay but with 10% more area compared to 
the solutions obtained by TDO, Furthermore, TDO spent much less time in achieving 
these high quahty solutions. 

A comparison was performed between TDO and the RERE algorithm (Pan, P., 
Performance-Driven Integration of Retiming and Resynthesis, in Proc. Design 



-31- 




Automation Conf,, pages 243-246, 1999) on the combinational circuits. The RERE 
algorithm performs retiming for sequential circuits. On average, the solutions generated 
by RERE have 19% more delay and 25% more area after the optimization, and have 6% 
more delay and 25% more area than the solutions obtained by TDO. 

It is concluded that the TDO method is able to achieve solutions with superior 
quaUties compared with the state-of-the-art timing optimization algorithms. 

A further interest is to understand the impact of the timing optimization on the 
circuit delay after the layout design. Although the timing optimization is very effective 
in optimizing the circuit depth during the technology independent stage, it may cause 
the circuit delay to increase if the layout effect is ignored during optimization. 
Therefore, the layout-driven synthesis flow was described in the Layout-Driven Timing 
Optimization section above. 

A comparison of the area optimization (optarea)^ timing optimization (TDO) 
and clustering-driven timing optimization (layoutDrivenTDO) was performed. 
Although there is no dramatic delay reduction, layoutDrivenTDO is able to achieve 
better area and better delay results compared to TDO which does not consider the layout 
effect. 

Finally, a comparison was performed between the synthesis flow 
(layoutDrivenTDO) and Quartus. On average, the synthesis flow obtains 16% better 
delay results and 10% better area results compared to the Quartus results. 

CONCLUSIONS AND DISCUSSIONS 
The novel methodology for the general timing-driven iterative refinement-based 
approach has the following characteristics and advantages: 

1. It reduces or eliminates the need for the (incremental) timing analysis 
during the iterative refinement procedure. 

2. It allows the local transformation to be able to see a more accurate timing 
information (e.g., the arrival times at the transitive fanin signals) so that 
the transformations can be conducted in a more meaningful way to 
reduce the circuit delay. 
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3. Its flexibility makes it much easier to consider the design preferences 
(e.g., in hybrid FPGAs with both LUT clusters and Pterm blocks, it is 
better to use the same logic resources consecutively on a critical path so 
that the subsequent clustering procedure can pack these implementations 

5 into one cluster to reduce the circuit delay.) 

4. It provides a general framework to integrate several types of local 
transformations, such as logic resynthesis, mapping, clustering, and so 
on, to enable an integration of currently separated design processes. 

The TDO method integrates the novel mechanism for the general iterative 
10 refinement flow and the area recovery technique using the restrictive iterative 
resubstitution. It is generally able to outperform the state-of-the-art algorithms 
consistently while significantly reducing the run time. 

Specific blocks, flows, devices, functions and modules may have been set forth. 
However, one of ordinary skill in the art will realize that there are many ways to partition 
15 the system of the present invention, and that there are many parts, components, flows, 

modules or functions that may be substituted for those listed above. 

MThile the above detailed description has shown, described, and pointed out the 
fundamental novel features of the invention as applied to various embodiments, it will be 
understood that various omissions and substitutions and changes in the form and details of 
20 the system illustrated may be made by those skilled in the art, without departing from the 
intent of the invention. 
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